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Abstract 

Dirac operator eigenvalues split into two when subjected to two different external vector sources. 
In a specific finite- volume scaling regime of gauge theories with fermions, this problem can be 
mapped to a chiral Random Two-Matrix Theory. We derive analytical expressions to leading order 
in the associated finite-volume expansion, showing how individual Dirac eigenvalue distributions 
and their correlations equivalently can be computed directly from the effective chiral Lagrangian 
in the epsilon-regime. Because of its equivalence to chiral Random Two-Matrix Theory, we use 
the latter for all explicit computations. On the mathematical side, we define and determine gap 
probabilities and individual eigenvalue distributions in that theory at finite TV, and also derive the 
relevant scaling limit as N is taken to infinity. In particular, the gap probability for one Dirac 
eigenvalue is given in terms of a new kernel that depends on the external vector source. This 
expression may give a new and simple way of determining the pion decay constant F w from lattice 
gauge theory simulations. 



1 Introduction 



One of the most challenging - and perhaps most interesting - problems associated with lattice gauge 
theory simulations of QCD is that of the chiral limit. Based on a variety of different approaches it is 
now possible to perform numerical simulations of the theory with two light dynamical quark flavours, 
at least in modest space-time volumes. By "light" quarks we mean quarks that are very close to the 
actual physical masses of the u and d quarks in QCD. Even if the masses of the physical u and d quarks 
had turned out to be much heavier (on the typical QCD scale Aqcd), one would like to explore the 
chiral limit of the theory in it own right. This is because the theory in this limit separates into to 
two disjoint regimes, of which the low-energy part can be treated in a systematic manner by means 
of effective field theory: the chiral Lagrangian based on the spontaneous breaking of chiral symmetry. 
This low-energy theory of QCD with very light quarks can be understood in much the same way that 
the low-energy limit of QCD without quarks matches on to an effective string theory description, and 
both limits are of interest. 

The so-called e-regime of QCD pQ is particularly useful for studying the chiral limit of QCD in 
finite volume. It is well known how a universality class of chiral Random Matrix Theory [2] provides 
an intriguing alternative description of the leading-order expressions for Dirac operator eigenvalue 
correlation functions in this regime, results that also can be derived directly from the low-energy 
effective field theory [3j dj. Even the distributions of individual Dirac operator eigenvalues follow 
from a systematic expansion in the chiral Lagrangian framework [5]. All of these analytical results 
depend on just one single low-energy constant of QCD, that of the infinite-volume chiral condensate 
S. From a lattice gauge theory viewpoint, this provides a new and unusual way of determining this 
low-energy constant of QCD by measuring the lowest-lying Dirac operator eigenvalues. For some 
numerical analyses see, e.g., refs. [El [7]. 

Recently, a new scheme was proposed which uses Dirac operator eigenvalues for determining the 
pion decay constant F n in a somewhat similar manner [8]. Based on the chiral Lagrangian formulation, 
the suggested method made use of a spectral 2-point function associated with two different Dirac 
operators, 

A# = [jp{A) + iw^ = a^4 n K (i.i) 

corresponding, in the case \i = \x\ = to imaginary isospin chemical potential. Equivalently, the 
two Dirac operators (II. lh are simply in a constant background Abelian gauge field, but with different 
"charges" . In the e-regime, the chemical potential (i couples directly to F n in the form of the finite- 
volume scaling variable (1 = fiF n y/V. Because the sensitivity to \x is quite drastic for the spectral 
2-point function, this provides a clean method for extracting F n . There is sensitivity to the parameter 
\x (and hence F n ) also in other observables in the e-regime [9j [10]. Alternatively, one may use a real 
chemical potential to determine F n |llj . with the same finite- volume scaling. The Dirac spectrum is 
complex in that case. 

The chiral Lagrangian approach of ref. [8] can, to leading order in the e-regime, also be re-cast in 
terms of Random Matrix Theory, this time a Random Two-Matrix Theory [12] . All eigenvalue density 
correlations are equivalent in the two theories j8[ @]. One loop corrections to both £ pQ and F v |13j 
have been computed in the e-expansion. To that order they simply amount to finite-volume corrections 
to the infinite- volume quantities S and F n ; the effective theory otherwise remains unchanged. It is 
of course important to know the size of these finite-size corrections if one wishes to determine £ and 
Fn from the eigenvalues of the Dirac operator by means of lattice gauge theory simulations at finite 
volume. An alternative method for extracting £ and F^ in the e-regime of QCD can be based on fits to 
vector and axial vector two-point correlations functions; also here finite-volume corrections are know 
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analytically at sectors of fixed gauge field topology [14] . See the recent review [15] for a summary of 
these different approaches. 

In the Random Matrix Theory formulation analytical computations are substantially simplified, 
and in ref. [12] all possible spectral density correlation functions associated with the two Dirac opera- 
tors were found analytically. This includes all spectral functions in both the quenched and unquenched 
theory, and even all spectral correlation functions associated with "partially quenched" spectral cor- 
relation functions, where there is no back-reaction of the chemical potential on the gauge field config- 
urations. This latter set of spectral correlation functions give the most fruitful way of extracting F v 
from lattice data since one can make use of ordinary configurations without chemical potential. Once 
all spectral correlation functions are known, one should in principle have all spectral data, and thus 
be able to reconstruct individual eigenvalue distributions as well. Indeed, this was precisely what was 
found in ref. [5] for the case without imaginary chemical potential. In the first part of this paper we 
show that this is also the case here0 It was also shown in [18] that in the quenched theory in the limit 
of large chemical potential all spectral and individual eigenvalue correlations factorise into quenched 
one-matrix theory quantities. 

In the second part of this paper we aim at determining in as closed form as possible the precise 
analytical expressions for individual eigenvalues distributions. Such expressions may turn out to be 
very useful alternatives for extracting F n from lattice simulations in those cases where only a very 
small number of "good" eigenvalues are available, and where it therefore may be difficult to construct 
the spectral 2-point correlation function with good statistics. Here we concentrate on the smallest 
eigenvalue, but provide the general framework for computing others. Based on the Random Two- 
Matrix Theory representation, we derive an explicit and quite compact representation for any finite N. 
Taking the scaling limit with N — > oo, this provides the sought-for analytical expression for the lowest 
Dirac operator eigenvalue distribution in the appropriate finite-volume scaling regime. Remarkably, 
the final expression is not much more involved than the one without external vector sources. 

From the point of view of mathematical physics, the resulting solution for the distribution of the 
smallest eigenvalue in the chiral Random Two-Matrix Theory is of interest in its own right. The 
solution cannot be mapped on to an analogous one-matrix theory, so the distribution is new and 
presumably corresponds to a new universality class that is parametrised by one real number 5. For 
this reason we include a rather detailed derivation, even though the resulting formula is all that is 
needed for the purpose of applications to lattice gauge theory simulations. 

2 Eigenvalue correlations in chiral Random Two-Matrix Theory 

We start by giving the theory we will solve for individual eigenvalue correlations, chiral perturbation 
theory in the epsilon regime with imaginary chemical potential 

Z { u Nf) =f dU ^ etU y e kVFZTr[U,B][Ul,B]+^VTr(MlU+MUt) _ ^ 
Ju(Nf) 

Here the matrix B =diag(/iiljv 1 ,/X2lAr 2 ) contains the two different chemical potentials, and A4 = 
diag(mi, . . . , rriN f ) is the quark mass matrix of the N± + N2 = Nf flavours. 

This theory and all its spectral density correlation functions are completely equivalent to the 
chiral Two-Matrix Theory with imaginary chemical potential that was introduced in ref. [12]. The 
equivalence for the two-point function follows from [8], for all higher density correlations it was proven 

1 A preliminary account of this was presented at a conference last year [16] . A comparison with Monte Carlo data 
from lattice gauge theory was presented at the same meeting [17] . 
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in [4] . It is defined as 

Ni N 2 



Z { u Nf) ~ J d$d* e -^(* t *+* t «) detfPa+m/a] fj det[P 2 + m /2 ] , (2.2) 
where is given by 



/1=1 /2=1 



The operator remains anti-Hermitian because the chemical potentials are imaginary. Both $ and ^> 
are complex rectangular matrices of size N x (N + v), where both N and v are integers. The index v 
corresponds to gauge field topology in the usual way. 

Referring to ref. |12| for details, we immediately write down the corresponding eigenvalue repre- 
sentation: 

f00 N / Ni N 2 

Z » = / II dx i dy i (x i y i y+ 1 ]J (xf + m 2 fl ) ]J (y 2 + m%) 

J ° i=l \ fl=l /2=1 

x A N ({x 2 })A N ({y 2 }) det [I v (2dNx iyj )) e~ N ^ c ^+ c ^ , (2.4) 

up to an irrelevant normalisation factor. Here the Xj's and y^'s are the real and non-negative entries 
in the diagonal matrices X and Y, defined by 

$i = $ + = UiXVl , 

$ 2 = $ + = U 2 YV^ . (2.5) 

Because of this redefinition the matrices <E>j become coupled in the exponent, leading to the above 
structure after integration out the unitary matrices Ui and V{. This leads to he following combinations 
of the two chemical potentials in eq. (|2.4p : 

a = (l + n 2 2 )/5 2 , c 2 = (l + rf)/6 2 , 
d = (1 + m/ u 2 )/<5 2 , 1 - r = d 2 /(cic 2 ) , 

5 = fj, 2 - m , (2.6) 

where the latter will become useful in section HJ 

The joint probability distribution function which is proportional to the integrand in eq. (|2.4f) is 
defined as 



1 N ( N\ N 2 

vi Nf \{x},{y};{mi},{m2}) = ]J (xaft)^ 1 ]J (s? + m 2 ;i ) J] (vl + ™% } 

Zu i=l \ fl=l /2=1 

x A N ({x 2 })A N ({y 2 }) det [I v {2dNx iyj )\ e -"£f (2.7) 



where An({x 2 }) = Ylf>ii x 'j ~ x i) ls the Vandermonde determinant. It is normalised to unity 



poo N 

1= / l[dx i dy i vi Nf) ({x},{y};{m 1 },{m 2 }) . (2i 
J° i=i 
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From the joint probability distribution we can define an (n, k)-density correlation function 

,2 »™ N N 



R k ,i{i x }k,{y}i) = ^ _ fyi( N _ ^ fl fl ^'(WiMiWil-W) 

(2.9) 



i=fc+l j=J+l 

JV! 2 1 



(N-k)\(N-l)\ z (»f) 

,oo N N 



poo iy 1 v 

/ II dxi II d % det A JV ({x 2 })A iV ({y 2 }), 



i=k+l j=l+l 

where we have moved the exponential and masses into the determinant, introducing 

N\ N 2 

wl N '\xi, yj ) = (x tyj Y +1 e- N ^ +c ^I u (2dNx iyj ) \\ (x\ + m%) J] (Vj+m%) . (2.10) 

/1=1 /2=1 

Obviously i?o,o = 1 is normalised to unity. The Rk,i({%}k> {y}l) can be expressed in terms of a 
determinant of four different kernels. These are given by the (bi-)orthogonal polynomials and their 
integral transforms with respect to the weight function eq. (|2,10p . and we refer to [12] for details. The 
Rk,i({x}k, {y}i) w dl be the building blocks to compute the gap probabilities as well as the distributions 
of individual eigenvalues of both type x and y. 
We define the following gap probabilities as 

jy|2 r-s roo r-t roo 

Ek,i(s,t) = (jy _ fcWjV - 11) J dxi --- dxk J dx k+1 ...dx N J dy 1 ...dyij dy l+1 ...dy N 

x vi Nf \{x},{y};{m 1 },{m 2 }) , for k, I = 0, 1, . . . , N , (2.11) 

where in the sequel we suppress the dependence on masses and topology for simplicity. The Ek,i(s,t) 
give the probability for general k, I G {0, . . . , N} that the interval [0, s] is occupied by k x-eigenvalues 
of V\ and [s, oo) is occupied by (N — k) x-eigenvalues, and that the interval [0, t] is occupied by I 
y-eigenvalues of T>2 and [t, oo) is occupied by (N — I) y-eigenvalues. It also depends on the masses and 
on |Ui 2 which we have suppressed here. 

Similarly we can define the probability to find the A;-th x-eigenvalue at value Xk = s, and the l-th 
y-eigenvalue at value y\ = t, to be 



Pk,i( s ^) = k \k^) l {j) J dxi ' ' ' dxk - 1 j g dx k+i---dx N J dy t . . . dyi-x dy l+1 ...dy N 

x vi Nf \xi,.. .,x k _ x ,x k = s,x k+1 , ...,x N ,yi,.. -,yi-i,yi = t,y i+1 , . . . , y N ; {mi}, {m 2 }). 

(2.12) 

Here the eigenvalues are ordered, x\ < . . . < xn and y\ < . . . < yN, and obviously k, I > 1. The fact 
that the p k} i(s,t) are probabilities that are normalised as 



ds / dtp k)l (s,t) = 1 , (2.13) 



can be seen along the same lines as for a single set of Dirac operator eigenvalues, as was shown in the 
appendix of [S]. 
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Because we have two sets of eigenvalues we may also define mixed gap-probability distributions as 
they will occur in intermediate steps. There are two different functions defined as 



Nl fN\ f s f°° [* f°° 

Ep k ,i(s,t) = 7]y _ k y l i l J J dx 1 ...dx k j dx k+1 ...dx N J dy 1 ...dyi_ 1 j dy i+l ...dy N 

x vi Nf \{x},yi, . . . ,yi-i,yi = t,y i+1 , . . . ,y N ;{m 1 },{m 2 }) , 

for k = 0, ...,N &ndl = l,...,N , (2.14) 



fN\ Nl f s f°° r f°° 

pE k ,i{s,t) = kl^j -—— J dxi...dx k -ij dx k+ i . . . dx N J dy x ...dyij dy t+1 ...dy N 

x 7>i Nf \xi,...,x k „i,x k = s,x k+1 ,...,x N ,{y};{m 1 },{m 2 }) , 

for k = 1,...,N and I = 0,...,N . (2.15) 

The first quantity eq. (|2.14p gives the probability that [0, s] is occupied by k of the x-eigenvalues of 
T>\ and [s, oo) is occupied by (N — k) of the x-eigenvalues, given that yi = t, where the y-eigenvalues 
are ordered. The second quantity eq. (|2. 15|) gives the probability that [0, t] is occupied by I of the 
y-eigenvalues of V 2 and [i, oo) is occupied by [N — I) of the y-eigenvalues, given that x k = s, where 
again the x-eigenvalues are ordered. 

These definitions include for example the probability pE k ^(s,t = 0) to find an eigenvalue of the 
first type at x = s, where all y-eigenvalues are integrated out. We will return to this in section |H 



3 Gap probabilities and individual eigenvalues from densities 

We use the simple identity 



{a -b)> = £(-!)' ] V'-^ , (3.1) 



and choose a = dx and b = dx to replace all the (N — k) cfa>integrals dx in eq. (|2.11|) by 
a — b, and likewise for the corresponding y- integrations. We obtain 



AH 2 f 
fc ' KM) " (N-k)\(N-l)l J 



dx t ... dx k / dyi . . . dy, 



n 







N-k /jy_iA / roo\ N-k-i / ,»> x ' 



N ~ l /AT — 7\ / poo\ N-l-j / r t\j 

E^yy j )(l ) [JJ dyi+i---dyNVr\{x},{y};{m 1 },{m 2 }) 

Yl — ttj — / dx x ...dx k+i / dyi...dyi +j R k +i,i+j(xi,...,x k+i , yi,...,y l+j ). 
7^ T^n * ! J ! Jo Jo 



i=0 j=0 



(3.2) 



Here we have used the invariance of the joint probability distribution under permutations of both 
{x} and {y}. The formula (|3.2p neatly expresses the gap probability in terms of spectral correlation 
functions of both sets of eigenvalues. The latter can be derived from /c-point resolvents in chiral 
perturbation theory by enlarging eq. (|2.ip to the corresponding supergroup integral, see [3|. Thus we 
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have shown how in this setting also gap probabilities and individual eigenvalue distributions follow at 
the level of the chiral Lagrangian. 

We may introduce a generating functional for all gap probabilities, 



E(s,t;Z, V ) = ^2(-^y(- v y— d Xl ...dxi / , 

„■ ~— n • JO Jo 



N 

i ; / 

dx\...dxj R id (x 1 ,...,x i , yi,...,yj), (3.3) 



i,j=0 

where the term at i = j = gives unity. It immediately follows that 

Qk Ql 



E kjl (s,t) = (-l) fc+ '^^T £-.../:i-„) 



, for fc,i = 0,1,..., iV . (3.4) 

?=l,r?=l 



We will now relate gap probabilities, mixed and individual eigenvalue distributions to density corre- 
lations. It can be easily shown that 

-^-E k;l (s,t) = k\(pE ktl (s,t)- P E k+ljl (s,t)) . (3.5) 

For k = I = we have 

—E , (s,t) = - P E h0 (s,t) , (3.6) 
as from the definition pE k i has k > 1, and thus we set pEoj(s,t) = 0. Similarly it follows 

^-E k)l (s,t) = l\(Ep k> i(s,t) - Ep k> i +1 (s,t)) , (3.7) 
at 

where again Ep k fi(s,t) = 0. If we differentiate the mixed correlators we obtain 

d 

-Q^Ep kj i(s,t) = kl(p kj i{s,t) -p k+lt i(s,t)) , (3.8) 
d 

Q^pE k ,i(s,t) = l\(p kt i(s,t) -Pk,i+x(s,t)) . (3.9) 
Finally, if we differentiate the gap probabilities twice we arrive at 

d 2 



dsdt 



E k i{s,t) = k\l\(p k As,t)-p k+ ii(s,t)-p k i + x(s,t)+p k+ ii + i(s,t)) . (3.10) 



Of course the order of differentiation does not matter, as one can easily convince oneself. The boundary 
conditions to be imposed here and in eqs. (|3.8p and (|3.9p follow from 



d 2 

;E k , (s,t) = - k\(p k>1 {s,t) -Pk+i,i(s,t)) 



dsdt 
8 2 

—E 0tl (s,t) = - H(pi,i(s,t)-p 1)l+1 (s,t)) . (3.11) 

Again from the definitions we have p k fi( s i~t) = Po,i( s ,t) = 0. 

The probabilities p k i(s,t) can be solved for the (mixed) gap probabilities in three different ways. 
Summing over k in eq. (|3.8p . or over I in eq. (|3.9j) we obtain 

Pn+l,l(s,t) = ~^2y^ E Pk,l{s,t) , 
n 1 d 

Pk,n+l{s,t) = J\'Q t pEk ^ S ^ ■ ( 3 - 12 ) 

1=0 
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Alternatively one can sum over both k and I in eq. A3. lQf) to obtain an expression in terms of gap 
probabilities alone 

n q 1 d 2 

P„+l,,+l(s,t) = +J2Y,^Vd^- t E kA S ^)- ( 3 - 13 ) 
k=0 1=0 

Let us give some examples. For the simplest case of k = I = we get the probability that the 
interval [0, s] is free of x-, and the interval [0, t] free of y-eigenvalues: 

/oo 
JJdxi 11 d Vj vi Nf \{x}, {y}; {mi}, {m 2 }) , (3.14) 
- i=i Jt j=i 

and we obtain 

d 2 

—E 0>0 (s,t) = p h i(s,t) . (3.15) 
Explicitly we have for this gap probability the expansion eq. (|3.2[) given already in [16] 

Eo.o{s,t) = 1- dxRi >0 (x)- / dyR Q> i(y) + - [ dx x dx 2 i?2,o(^i, x 2 ) 
Jo Jo 1 Jo 

1 /■* 

+ o / dyidy 2 Ro, 2 {yi,y 2 ) + ... 
z Jo 

+ dx dyRi t i(x,y) - - / dx\dx 2 I dy R 2 ,i(xi, x 2 , y) 
Jo Jo 1 Jo Jo 

~o I dx I ^1^2-^1,2(2,2/1,2/2) + ■•• ■ (3.16) 
1 Jo Jo 

The terms in the first two lines containing only s- or i-dependent integrals are annihilated by the 
differentiation in eq. (I3.15p . and we get to the same order 



pi.i(s,t) = Ri t i(s,t) - / dx R 2j i(x,s,t) - / dy Ri j2 (s,t,y) + ... . 

Jo Jo 

4 An exact expression for the first eigenvalue distribution 



(3.17) 



In this section we derive a closed expression for an individual eigenvalue distribution, the probability 
to find the first eigenvalue of T>\ at s irrespective of the position of the ^-eigenvalues. Our solution 
given in terms of a new kernel and polynomials holds for any number of flavours Ni and ^2 and 
arbitrary chemical potentials fM\ and \i 2 . In particular, we can partially quench the type-1 flavours 
(putting N\ = 0) with fi\ 7^ in gauge theory with type-2 physical sea-quark flavours (i.e.,N 2 ^ 0) 
with fx 2 = 0. This case is probably the most interesting for applications to lattice QCD. 



4.1 The finite- iV solution 

We first consider the gap probability that the interval [0, s] is empty of x-eigenvalues, 

/oo poo . . 

dx x . . . dx N / dyi... dy N VI f) ({x}, {y}; {mi}, {m 2 }) . (4.1) 
JO 

From this the sought probability follows by differentiation, pEi t o(s,t = 0) = —d s Eofi(s,t = 0). We 
will now perform a series of steps before arriving at an exact expression for finite N. The appropriate 
large- N scaling limit will be taken in the next subsection. 
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First, recalling the definition of the joint probability distribution eq. (12. 7p we can use the fact 
that all y-eigenvalues are integrated out in eq. (|4.ip for a symmetry argument. The Vandermonde 
determinant Aj^({y 2 }) and the determinant of the Bessel function are antisymmetric. Therefore we 
can replace the latter by its diagonal part times N\ 



£o,oM) 



dx 1 ...dx N / dyi . . . dy N A N ({x 2 })A N ({y 2 }) 



m 

N / Ni N 2 

II i*iVi)" +1 e- N ^ +c '^I u (2dNx iyt ) [J (x? [] (yf + 

i=l I /1=1 /2=1 



"1/2) 



(4.2) 



In the next step we use a known identity relating the Laguerre weight times the /-Bessel function to 
an infinite sum over Laguerre polynomials (see e.g. eq. (B.7) in [12]). With this decomposition we 
can exploit the orthogonality properties of these polynomials. For simplicity of the proof we will first 
consider JVi = N 2 = 1 (Nf = 2) with masses mi and 7712, and later give the general result for any 
numbers of flavours. We thus have 



£o,oM) 



7(1+1) 



N 



dxi... dx N / dyi... dy N A N ({x 2 })A N ({y 2 }) JJ(xf + m\){y 2 + m|) 



8=1 



A" 



X 



n 



i=l 



r ~ { x iyi ) 2 » +1 e- NT ^ +c ^Y. 

n t =0 



nj\(l - tY* 



L» ni {NT Cl x 2 )L» ni (NTC 2 y 2 )y 

(4.3) 



with r = 1 — d 2 /(C1C2). Next we include the mass 777,2 into the Vandermonde determinant A]\f({y 2 }), 
and then replace it by a determinant of Laguerre polynomials normalised to be monic 

N 

A N ({y 2 })H(y 2 + m 2 ) = A N+1 ((im 2 ) 2 ,{y 2 }) = det [(-^(A^r^W)] , (4.4) 

3,k=0,...,N 

i=l 

where yo = im%. We observe that the Laguerre polynomials L^{Nrc 2 yk) now all appear with their 
corresponding weight function y 2u+1 e~ NrC2V ^ , except for yo of cours^l. Writing the determinant eq. 
(|4.4p as a sum over permutations we can integrate out all variables y\ to yjv successively, each integral 
killing one infinite sum over n,. Thus, under the permutation each Ln} \NTC2y 2 ) gets replaced by 
Ln} {Ntc\x 2 ) times the norm from the integration and the remaining factor from inside the sum. 
We can therefore rewrite the result again as a determinant, with the first row containing the mass 
yo = im2 unchanged: 



£o,oM) 



{Nd) Nv T N { v+l) / ^ dxNAN{{x 2 }) JJ (x 2 + m?) x 2u+l 
J s „-_i 

Ll{Nrc 2 {im 2 ) 2 ) • • • ^^L^NT^im^ 2 ) 



-(1+1) 



-Ntc\x\ 



' L^rc.x 2 ) ••• MS^ L< f\Nrc lX 2 ) 



2(Ntc 2 )»+ 1 
2(Nrl 2 )»+ i Lq(NtCiX 2 n 



2(NTC 2 y 
2(Nrc 2 y 



HiMk^ L y {NTClx 2 



NJ 



(~) N Nl TV 



iNrC2 )N L N( NTC 2( lm 2 
-) N NU 

2(Nrc 2 ) N + v + 1 



(-rm^rr L (») {NTcixl) 



2(Ntc 2 ) n + v + 1 N l ivTL l x AT/ 



2 Note the additional r in the exponent comparing eq. (|4.2p and (|4.3p . coming from the identity for the /-Bessel 
function. 
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N\{Nd) N "T N ("+V nf =0 (l - Ty(NTc 2 ) 



(1+1) 



2 n (Ntc 2 ) n ^+ 1 ) 



N 



dx! . . . dx N A N ({x 2 }) J](xf + m\) xf +1 e- NTC ^i 



LUNTc 2 (im 2 ) 2 ) 
L v a {Nrc x xl) 

L»(Nt Ci x 2 n ) 



T)~^L^{Nrc 2 {im 2 ) 2 ) 
L^Nrax 2 ) 

mNr Cl x 2 N ) 



(1 



i=l 

-N TV 



t)- n L v N {Nrc 2 {im 2 ) 2 ) 
L» n (Nt Ci x 2 ) 

L» n (Nt Ci x 2 n ) 



(4.5) 



taking out common factors. Here we have denned the following notation for monic Laguerre polyno- 
mials 

L-(x) = (-l)"n! L v Jx) = > Y-ir* •••.\ , :/.; J L,, ^ • (4-6) 



Vf nn+j n!(n + t/)! ^ 
' (n-3)\(y + j)\j\ 



The last determinant in eq. (j4.5h can be almost mapped to a Vandermonde determinant, using the 
following identity: 



Lg(Mf) • 




L» N (X 2 ) 




•• L»(Xl) ■ 


-ko(^v) " 


■■ L){X%) ■ 


l n( x n) 


1 


x 2j 


T N TV ( M l\ 

X 2N 


1 


N 


v2N 



where we have defined 

M 2 = NTc 2 {im 2 ) 2 and X\ = Ntcix\ . 
A proof of this relation is given in appendix lAl 



We can now change variables Xj 



u i 



x 2 , and perform the shift Uj 



u i 



(4.7) 



(4.8) 



s to obtain 



integrations dzj in eq. (|4.5p . The change of variables and subsequent shifts induce the following 
changes: 



A N ({u}) 
(uj + m\) 



U ; 



-NrciUi 



Aiv({4) , 

(zj + s 2 + m\) = 

( Zj +s 2 r , 

e -NTCiZi e -Nrcis 2 



{ Zj + m[ 2 ) 



(4.9) 



In other words: the Vandermonde determinant remains invariant, the mass rrt\ receives a shift to 
s 2 + m 2 , the topology term becomes a i^-fold degenerate mass term with mass s 2 , and the 
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weight is shifted by a constant factor. While this is just as in the chiral one-matrix theory, the 
difference here is that the almost Vandermonde eq. (|4.7p is not invariant, and becomes 





r J jv I M$ nj 


•• (i_ T )iv^M r 


1 


Af 


X 2N 


1 


N 


v2N 
N 



(4.10) 
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{NT Cl {z 1 + s 2 )y 



(nt CI ( zn + s 2 )y 



T N TV ( M l\ 

(1_ T )IV-^M T ) 



(Nt Ci (z n + s 2 )) N 

T TV( 1V1 2 \( C2 

4=0 (l-rf^l { — ){ D 



Z N 



zf 



yN 

N 



where we have defined 

S 2 = Ntcis 2 and Z k = Nrciz k . (4.11) 

Here we have subsequently added column^) to turn back to monic powers in Z k . The fact that it is 
not quite invariant illustrates the fundamental property that there is apparently no way to map the 
present two-matrix problem onto an equivalent one-matrix problem. 

In the next step we turn the monic powers Z? back into monic Laguerre polynomials, using again 
the invariance of the determinant. This will introduce yet another sum over the Laguerre polynomials 
in the mass variable M| in the first row. Because of the shift eq. (|4.9p we only need to generate 
Laguerre polynomials with topology v = 0, Lj =0 = Lj, to obtain the polynomials orthogonal to the 
shifted weight. We thus obtain for the last determinant in eq. (|4.10p 



Qo(Mf) 
L (Zi) 

Lq(Zn) 



Qj(Ml 

Hz,) 



Qn(M 2 ) 
L N (Z 1 ) 

Ln(Zn) 



(4.12) 



which defines new polynomials 



Qn{M 2 ) = 



^ (n 



(_)n+i( n J)2 



1=0 v / \ / 

n I 



V r 



-) k+i i\ 



-S 2 ) j 



(~) n n 



=o 

n 

fc=0 v 



k=0 



(l-k)\k\(l 



:L k (M 2 



- k L k (M 2 )L-\{-S 2 



(4.13) 



In the first step we have swapped sums, ^j=o Yll=o ~^ J2?=o Sj=z> sucn that the powers in S 2 give 
the Laguerre polynomial L l J— S 2 ). In the second step we have used the identity (|A,8P backwards in 
order to take the argument 1/r out of the first Laguerre polynomial, in choosing z = M| and w = 1/r 
in eq. (|A.8p . This goes at the expense of introducing another sum. After swapping again sums to 
Sfc=o Sr=fc> t ne latter sum over the generalised Laguerre polynomial in — S 2 can be simplified, using 
the following identity, 

n—k , . ,\ 



l: 



-s 2 



=o 



(4.14) 



3 Usually one goes from monic powers to monic polynomials in this step. It is easy to invert this by defining Z[ 
and then going from monic powers Z'^ to polynomials (Z[ — S ) J = Z\. 



Z,+S 2 
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A proof of this simple identity is presented in the appendix[B] Note that all polynomials L n ^ k (—S 2 ) 
are proportional to — S 2 , except for n = k as Lg^-S 2 ) = 1. This will become important when 
computing the normalisation in the limit s — > 0. 

The explicit appearance of a new set of polynomials is again a reminder that we cannot map the 
problem onto a one-matrix problem. Of course, in the limit of the deformation parameters fix 2 —* 0, 
the polynomials Q n (M 2 ) reduce to Laguerre polynomials. In the form given in the last line of eq. 
(14.131) the new polynomials Q n are amenable to the large-iV scaling limit that we take in the next 
subsection. 

Returning to the gap probability, in the last step we now replace in eq. (14. 5h the determinant 
Aat({x 2 }) = A]\r({z}) times the mass term by a larger Vandermonde determinant, and then replace 
monic powers by Laguerre polynomials monic in the arguments Zk- The v degenerate masses obtained 
after the shift eq. (|4.9p can be dealt with by first taking them different, and then taking limits by 
multiple application of l'Hopital's rule. For simplicity we set v = in all of the following. We have 



N 



Aat({4)I[( 



Zi + m' 2 ) 



AAr+i^imi) 2 ,-^}) 



i=l 



det [(-yjliNrci^LjiNrdZk)] , (4.15) 

J,K — U,...,iV 



where we define zq = imj. 

Let us collect what we have derived so far: 



£o,oM) 



L (M( 2 ) 
Ln(Z,) 



^!nf=o( 1 - r)>{NTcz)-i{NTcf t 



-N 2 TC 1S 2 



r(l+l) 



/2\ 



2 2 ^(A^c 2 ; 



A' 



, dzN e 



Lj(M[ 



L N (M{ 2 ) 
Ln(Zx) 



L (Z 



N 



L 3 {Z 



Here the definition 



AT 



M[ 2 = NTCxiim'x) 



L N (Z 



N 1 



1 \2 



Qo(M|) 
L (Zx) 

Lq{Zn) 



-Nrcximl + s 2 



Qj(Ml) 
Lj(Zx) 

Lj(Zjsi) 



Qn(M 2 ) 
L n {Zx) 



L N (Z 



N , 



(4.16) 



(4.17) 



has been used. We can now apply the orthogonality of the Laguerre polynomials with respect to the 
weight e ~ NTC ^ z i to compute the iV-fold integral over the determinants, applying the standard Dyson 
Theorem. We thus obtain the final answer for finite N: 



£o,oM) 
K N+1 (M[ 2 ,M 2 ) 



Ce 

A' 



-N 2 



TC\S 



K N+1 (M[ 2 ,M 2 



Y^ { -fHM[ 2 )Q,(M 2 ) 



(4.18) 



j=0 



This result defines a new kernel of the polynomials Lj and Qj in the (shifted) masses. The constant C 

i/E 



A 



1 



r(l+l) 



jj=Q (x-tV iMx)Lj (M2 ) that is inversely proportional to the partition function Zq ' ensures 
the correct normalisation lim s _,o Eq ; o(s, 0) = 1. It can be obtained independently by computing the 
partition function 



r(i+l) 



nf = o(j') 2 (l - Ty(NT Cl )-l(NTC 2 ) 

2 2N (Ncx) n (Ntc 2 ) n 



A 



3=0 



(l-r> 



-l^mDl^mI 



(4.19) 
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following the same steps as before but setting s = 0. The calculation simplifies in eq. (I4,10p so that 
the Qj's become Laguerre polynomials. Indeed as a check we can take 



lim Q n {M 2 ) 



:2 



(-) n n! 



(4.20) 



1 contributes to the sum. This already indicates that this last 



where only the term L~^_ k=Q (— S 2 
term in the sum is special. 

As a further independent check we may take the limit fj,i, )X2 — > (=> r — ► 0). In this limit the two 
Dirac operators become equal, T>\ = T>2- and we should recover the known one- matrix theory result 
[20] . Indeed, we get 



lim — - 

Mi>A»a->0 n\ 



-Qn{Ml 



^TL k (-Nm 2 2 )L-\(-Ns 2 ) 

k=0 



(4.21) 



^ (n -?')!(/' ) 2 ^ 

j=0 y ■>> u 1 l=Q 



L, 



N{ml + s 1 



i.e., Laguerre polynomials of shifted mass just as for the first flavour m! x . This follows from the first 
and last line of the definition eq. (|4.13p . Inserted into the kernel eq. (j4.18j) we obtain the one- matrix 
theory result for the gap probability in terms of the partition function of 2 flavours with shifted masses. 

It is straightforward to see that for more flavours, N\ > 1 and N<i > 1, the very same steps still 
go through (see also the corresponding determinant identity in appendix |Aj). The only difference is 
that the there will be more rows with masses of flavour N\ in the determinant eq. (14. 15ft . and more 
masses of flavour Ni in the new polynomials Qj in eq. (|4.12p . The absorption of the mass terms into 
a larger Vandermonde determinant leads to inverse Vandermonde determinants in each of the N\ and 
N2 masses, which can be taken out of the integral. We arrive at 



£o,oM) 



-N 2 tc\s- 



POO 

/ dz x ...dz N e~ N ^=° TClZi 
Jo 



(4.22) 



U{M' 2 



Lo(Zi) 



LN+Nl~l(Mj: 1=l ) 

ijV+ATi-l^Jv) 



Qo(M 



/2=1/ 



Qo(m 2 N2 : 



Qn+n 2 -i(M 2 2=1 

Qn+n 2 -i{M 2 N2 ) 
Ln+n 2 ~i(Zi) 



L 



N+N2-1 



(Z. 



Nj 



The orthogonality of Laguerre polynomials can be exploited in the manner of ref. [19]. This leads to 
the following determinant expressions. For an equal number of flavours N% = N2 we hav^l 



Eq,o(s, 0) = const. 



-N 2 



f ) A, 1 ({ ra y ) A, 2 ({ m y) 



det [K N+Nl -i(M' f lM%)] 
1</1,/2<JVi 1 J 



(4.23) 

In the case where N\{N2) is larger, the determinant is of size Ni(N%) and contains additional polyno- 
mials Lj(M'jf) (Qj(M% 2 )) to fill up the additional columns (rows) [19] . 



4 We omit all mass dependent normalisation constants here that can be obtained easily. In particular they will cancel 
the Vandermonde determinants of the masses, see e.g. |12] for the partition functions. 



13 



Nx > N 2 : E 0>0 (s, 0) = const.^, e -ivW 2 ( 424 ) 

Z ( Nf) A Nl ({m 2 fl })A N2 {{m} 2 }) 

x det [Kn+N^M'jl, M| 2=1 ) • • • K N+N ^{M' f l, M 2 N2 ) L N+N , 2 (M' } 2 ) ■ ■ ■ L N+Nl ^{M' f \)\ 
where we display the /1-th row in mass M' f l of flavour Ni , and 

N 2 > Ni : E 0fi (s, 0) = const.— — g -ivW 2 (4 25) 

x det [K N+Nl ^(M] 2 , M' f 2 1=l ) ■ ■ ■ Kn+n^M^M'^) Q N+Nl {M] 2 ) • • • Q N+N2 ^(M] 2 )] 

Here we have transposed the matrix to display the /2-th row in mass M% of flavour N 2 . 

For example, this includes in particular the interesting case of quenching the first flavour (i.e. 
putting Ni = 0) while keeping its chemical potential nonzero, [i\ 7^ 0. This quenched flavour can then 
be measured in the background oi N 2 = 2 flavours with masses mi and m 2 , 

e _ N 2 TCls 2 



N 2 = 2, N\ = : £ ,o (s, 0) = const— ^ det 

Z\ >(ml-m{) 



Q N (-NTc 2 mf) Q N+ i{-NTc 2 m\) 
Q N (-NTc 2 ml) Q N+ i(-NTc 2 ml) 

(4.26) 

In particular, setting the chemical potential of the dynamical flavours N 2 to zero, u 2 = 0, will not 
eliminate the other chemical potential Hi 7^ (see eq. (j2.6|) ) . or reduce to a known one-matrix quantity. 
This non-trivial /i-dependence due to the valence quarks can serve as a clean way to measure the pion 
decay constant F n from gauge field ensembles generated with dynamical light quarks that carry no 
chemical potential. 

Finally, as was pointed out earlier, the probability corresponding to non-vanishing gauge field 
topology v 7^ can be introduced by adding v extra masses of T>i, and then taking them to be 
degenerate with value s 2 . 

4.2 The large-iV limit 

In this subsection we take the large- N scaling limit by the same rescaling as in ref. [12], to which we 
refer for more details. We first derive the limits of all building blocks needed for the general case, and 
then specify the fully explicit result in three examples in subsequent subsections. 

All eigenvalues, the gap and the masses are rescaled in the same way (as would be y = 2Ny), the 
usual microscopic limit 

x = 2Nx , s = 2Ns , 
rhf = 2Nm f , fi f = 2Nu 2 f for f = 1,2 , 

5 = fa -fa. (4-27) 

All scalings including the chemical potential keeping Nfi 2 fixed can be read off from the chiral La- 
grangian eq. (|2.ip . For the various constants containing the u j this implies the following scaling: 



lim (1 - r) = exp 



-rt5 2 
2 



where t = j/N, r = k/j , 



lim ret = 1 for / = 1,2 . (4.28) 

N— too 
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For the Laguerre polynomials the following scaling holds: 

lim Lj(M{ 2 = -NTCim?) = I {Vt fh[) , (4.29) 

N,j—rOO 



lim 



recalling m 2 = rh\ + s 2 . Special care has to be taken in the asymptotic of the new polynomial, 

H -Q,-(Mf = -Nrc 2 m 2 2 ) = .lkn ^ ^-1_ L ^ M 2 ) L^{-S 2 ) + j^-^L^M 2 ) -l\ 

+ e^ tS2 I (Vi m 2 ) . (4.30) 

Here we have to split off the s-independent part Lq 1 = 1, which is the single term surviving in the 
limit s — » 0, and which hence ensures the normalisability of the probability in that limit. Usually 
neglecting a single term when replacing a sum by an integral amounts to removing a quantity of 
measure zero, which should be irrelevant. However, in our case this is not true as the convergence in s 
is not uniform. Therefore we have to treat that term separately and find the "anomalous" Jo-term in 
the scaling limit. That this procedure is correct is checked by computing the normalisation before and 
after taking the large- N limit. This curious phenomenon together with the appearance of Laguerre 
polynomials Lj_ k leads to our new microscopic kernel. 

The final answer for the microscopic limit of the new kernel in eq. (|4.18p thus reads 

K s {m' x ,m 2 ) = lim ±-K N+1 (M[ 2 ,M 2 ) (4.31) 

N^oo iV 

= 2 f dTT 2 I {Tm\) [ dRR 1 = e^ R2T2p I Q (RTm 2 ) s h(sT y/l - R 2 ) 
Jo Jo y/l - R 2 

+ 2 f dTT ei T2 ^/ (Tm , 1 )/o(rm 2 ) , 
J o 

where we have changed to squared variables. This kernel can no longer be related to a single partition 
functions of shifted masses, as it was the case in the one-matrix theory |20| . 

Likewise we obtain for the normalisation constant which is proportional to the partition function, 



N i rl , 

J (Tmi)/o(rm 2 ) • (4.32) 



lim — V\ — 1 -—LJMhLAMh = 2 I dTTexp 
iV-oo N (1 - t)J 3K l) 3K 2) J 



-T 2 l 2 
2 



Note that the first mass rh\ is not shifted here, in contrast to the previous equation. Partition functions 
of more flavours follow easily given the building blocks above, together with the general expressions 
given in |19^ [T2] . We now have all ingredients to obtain all gap probabilities with any flavour content 
by inserting the asymptotic kernel eq. (|4.3ip and asymptotic polynomials eqs. (|4,29p and (|4.30p into 
the respective eqs. (|4.23p - (|4.25p . normalised by the corresponding partition function. 

In the following we give three simple examples that illustrate these very general expressions. In 
order to guide the eye we mostly display the distribution of the first eigenvalues versus the correspond- 
ing eigenvalue density it has to follow. For comparison to Lattice results the gap probability that we 
give explicitly may be even more useful as it allows for a binning independent comparison with data. 
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4.3 Two light flavours 

Let us first consider the gap probability corresponding to two flavours N\ = N 2 = 1, as given in eq. 
(|4,18[) . Collecting the formulae from above we obtain 

-1 . ... X" 1 



Jim^o.oM) = 4^,0 (».°) = [l dTTe^ 2 / (Tmi)/ (Tm 2 )j exp ^ 

x dTI^IoiTmfy) J dRR-^=L= e^ R2T2p I (RTm 2 ) s hisT^/l^^) 

+ j\TTe^ TH2 I Q {Tm[)I Q {Trh 2 )^ . (4.33) 

Once more we can perform an analytic check by taking 5 — > in order to go back to the known 
one-matrix quantity |20j. Using the following so-called Sonine integral identity [21] (that also follows 
from the large- N limit of identity eq. (I4,2ip ) 



1 x 



dx I (mx)h(s\/L - x 2 ) + I {m) = I ( \/ m 2 + s 2 ) , (4.34) 



x 



we obtain the known gap probability [20J as a ratio of a two-flavour partition function with shifted 
masses over one with unshifted onesj As an illustration, we show the distribution of the first eigen- 
values p s \ Q 1+1 (s, 0) = —dgE s q g 1+1 (s, 0) for different values of 5 and compare it to the corresponding 
densities in fig. [TJ From [12] we have for the density 

(l+l),^ _ a m Jo dttJoit^Ioitm,) dtt e^ 2&2 J (tx)I (tm 2 ) 

Pi,o \ x > - P\mm\ x ) x ~n 1*2^2 : , ; : : ; : ; > ^-^J 

fodtte** 5 Io(^i)Jo(im 2 ) 
where we have introduced the one- matrix model quenched density p® MM (x) from eq. (|4.37p below. 

0.5r 



0.4 




Figure 1: The eigenvalue density and first eigenvalue for Nf = 1 + 1 with imaginary chemical potential 
6=1 (low red), 3( middle green), and 10 (upper blue curve), at fixed quark masses rh\ = 3, m 2 = 4. 

At 6 = eq. (|4.35p coincides with the corresponding one-matrix model density eq. (|4.38p below. 
For 6 = 1 the curve is still close to this density, compare to fig. [2j For 5 1 the curves approach 

5 The remaining integral is elementary and gives a 2 x 2 determinant of Bessel functions. 
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the one-matrix quantities of one flavour Ni = 1 with mass rh\ = 3 (the flavour corresponding to the 
^-eigenvalues gets quenched), compare again to fig. [2] below. This fact can be seen analytically, by 
taking the limit 5 — > oo in eq. (|4.35p and doing a saddle point approximation, 

lim f dtt e\ t2p J (tx)I (tm 2 ) ~ <T 2 e^ 2 J (x)/o(r7i 2 ) • (4.36) 

S— »oo JO 

The (5-dependent integrals get replaced by their values at the upper limit t = 1, and we have also 
computed the subleading coefficient for later convenience. After cancelling common factors in eq. 
(|4.35p we obtain, as we should, the one-matrix density for one flavour with mass rh\ as in eq. (|4.40p 
below. We have checked that the same limit applies to the first eigenvalue distribution. 

0.5r 




Figure 2: The eigenvalue density and first eigenvalue of the one-matrix theory: two flavours with 
rh\ = 3, m2 = 4 (low red), one flavour with rh\ = 3 (middle blue), and the quenched case (upper 
black curve). 

For the comparison above we give the following known one-matrix quantities [2CH [22] that are 
displayed in figure The quenched density and its first eigenvalues read 

P?mmW = l{Jo{xf + Ji(xf) , V Q lMM (s) = \se"^ . (4.37) 

The massive two-flavour density is given by 

n (2) f *\ n Q a Jq 1 dttJ^hjtrm) J 1 dttJ (tx)Io(tm 2 ) 

P\mm\ x ) — P\mm\ x ) x r i 7777771 \ T - ~\ Kft.oo) 

J dttIo{tm 1 )Io{tm 2 ) 

as well as its first eigenvalue distribution by 

Here primed masses are shifted according to m i 2 = mf + s 2 . We also need the one-flavour density and 
its first eigenvalue 

/£!«(*) = Am(*) - **(*) J.'^y'*') , (4.40) 

lo(mi) 

(i) /~\ 1 - -is 2 h{m 2 ) 

P\mm(s) = ~ s e 4 JTTfTY ■ ( 4 - 41 ) 
I lo\mi) 
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4.4 Partial quenching 



0.5 




Figure 3: The eigenvalue density and first eigenvalue for Nf = + 1 with imaginary chemical potential 
5=1 (low red), 3( middle green), and 10 (upper blue curve) and fixed quark mass rh\ = 3. 



As another example we can consider the partially quenched gap probability eq. (|4.26p with N\ = 
and one single {N% = 1) flavour of mass mi, 



4°o,o(*' ) = ex P 



1 ~2 1 52 

--r — <r 

4 2 



Qs{rhi;t = 1) 
I (mi) 



(4.42) 



where the extra ^-dependent factor e ^ 2 / 2 comes from the partition function Z^ 1 ^ that normalises 
this gap probability Its derivative is shown in fig. [3] together with the corresponding density [12J 



(o+i) /-n Q (~\ 
Pifl i x ) = Pimm( x ) 



cxp 



. Jq{x) 
'J (mi) 



dTTe^ TH2 1 (Trhi)J Q {Tx) . 



(4.43) 



Again we recover the one- matrix density eq. (|4.40|) when setting 5 = in eq. (|4.43p . For 5 = 1 the 
curve is still close to this one-flavour one-matrix result, see fig. EJ and for 5 3> 1 the curves approach 
the quenched one-matrix density. This can again be checked analytically by taking 5 — > oo of eq. 
(I4.43|) and using eq. (|4.36|> . The exponentials cancel but the prefactor 1/5 2 makes the second term in 
eq. (|4.43p vanish, leading to the quenched result. Again these limits 5^0 and 5 — > oo can also be 
checked for the gap probability. 

As the last and probably most physically relevant example we consider the partially quenched case 
of Ni = and N2 = 2 flavours with two possibly non-degenerate masses mi and mi 



4°S 1 (*»°) = exp 



s 2 





Q s (rhi;t = 1) 


dtQs(ihr,t) 




det 




t=i 




Qs(rh 2 ;t = 1) 


dtQs{rh 2 ]t) 


t=i - 



J (mi)m 2 Ji(m2) - J (m 2 )miJi(mi) 



(4.44) 



When setting fa = to have the sea quarks of flavour A2 free of chemical potential we simply have 
5 = — jj-i 7^ 0. The limit of equal masses can also be taken at the expense of a further derivative within 
each determinant. As stressed before, this case should be particularly useful for lattice gauge theory 
simulations, since it corresponds to ordinary configurations without chemical potential. 
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0.5 




Figure 4: The eigenvalue density and first eigenvalue for Nf = + 2 with imaginary chemical potential 
5=1 (low red), 3 (middle green), and 10 (upper blue curve), at fixed quark masses of flavour 7V 2 
rh\ = 3, m 2 = 4. 



The comparison to the spectral density given bjH 



Pimm( £ ) ~ ex P 
U 2 t 2 



■If 



(mili(mi)l (m 2 ) - m 2 /o(?™i)Ii(m 2 )) 



(4.45) 



o 



dttez Jo(xt)Io(rhit)\ — Jo(m 2 )(xJi(x) + 5 Jq(x)) — m 2 /i(m 2 ) Jo(x] 



!- J dtte^ 2 * 2 J (xt)Io(m 2 t)^o(?"i)(£«/i(x) + 5 2 J (x)) + mih(rhi)J {a 

is shown in figlH The one-matrix model result with two flavours is again recovered by setting 5 = 
given by eq. (|4.38[> . For 5 = 1 the curve is close to the two- flavour one-matrix result, see fig. EJ and 
for i > 1 the curves approach the quenched one-matrix quantities, see fig. [2J This matching can 
once more be checked analytically by taking 5 — > oo of eq. (|4.45|) . Using eq. (|4.36|> as well as the 



. 8 



. 6 



0.4 



0.2 




Figure 5: The integrated lowest eigenvalue distribution 1 — Eq^(s) for 1 + 1 flavours eq. (14, 33ft (lower 
red) and + 2 flavours eq. (|4.44j) (upper black curve), both at masses rh\ = 3, m 2 = 4 and 5 = 3. 

6 There is a term missing in the bottom right of the 3x3 matrix in the relevant formula eq. (3.53) of [12] . 
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cancellation of the two terms proportional to 5 2 in the last two lines of eq. (|4.45|> leads again to a 
complete quenching of all flavours. The same can be checked for the gap probability. 

Finally we can also compare directly the gap probabilities in our two-matrix theory for 1 + 1 
flavours and + 2 partially quenched flavours, where in fig. [5] we show 1 — Eq^s). For 5 = 1 the 
difference is still small but it grows rapidly with increasing 5 since both theories converge towards 
different limits for 5 — > oo as was pointed out earlier. Equivalently this results into the following 
comparison for the densities and first eigenvalues shown in fig. [6l It should be noted here that the 

0.5r 




Figure 6: The density and first eigenvalue for 1 + 1 flavours (right blue) vs. + 2 flavours (left black 
curve), both at masses rh\ = 3, = 4 and fixed 5 = 3. 

quantity 1 — Eo t o(s) is the integrated lowest eigenvalue distribution, by some considered a convenient 
quantity for comparison with the lattice gauge theory data. 

We end this section by pointing out that there is no analogous computation of a //-dependent Dirac 
eigenvalue distribution in the fully quenched case. The mixed two-point spectral correlation function 
has non-trivial /i-dependence [8] (and this dependence allows for the determination of a quenched 
value of F n using this technique). But the one-point function is //-independent in chiral perturbation 
theory because it is generated by the addition of just one valence quark; it is //-independent to all 
orders in chiral perturbation theory because the valence pions do not carry net baryon charge. 



5 Conclusions and outlook 

The two main results of this paper are the following. We have shown how individual distributions of 
the lowest-lying eigenvalues of Dirac operators that are subjected to two different external Abelian 
vector potentials ("imaginary chemical potential") can be derived from field theory. The results have 
been given in terms of generalised gap probabilities from which the distributions can all be derived. 

To compute the gap probabilities from field theory one needs to know spectral correlation functions, 
all of which can be given a well-defined meaning in the field theoretical setting. In particular, in 
the scaling region known as the e-regime, these eigenvalue distributions can be derived from the 
corresponding effective theory, the chiral Lagrangian. To make these computations concrete, we have 
used the equivalent Random Two-Matrix Theory to derive the distribution of the lowest Dirac operator 
eigenvalue in the e-regime of QCD with imaginary chemical potential. As stressed in the introduction, 
these analytical formulas may provide a very convenient way of determining simultaneously the infinite- 
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volume chiral condensate £ and the pion decay constant F n by means of numerical simulations in lattice 
gauge theory. 

We have given explicit formulas for the lowest individual distribution in terms of a new kernel, 
both in the case of full QCD with imaginary chemical potential, and for the analogue of partially 
quenched QCD in which quarks are dynamical, but do not carry chemical potential. Especially the 
latter may provide the most useful formulation in terms of comparisons with numerical lattice data. 
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A A determinant identity 

In this appendix we prove the following identity for any number of N2 flavours 
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Here we use the notation L v n (x) = x n + . . . for the Laguerre polynomials in monic normalisation 

K(x) = (-l)«n! K (x) = E(-l)^ "^L, * 3 ■ ( A - 2 ) 



3=0 



For simplicity we will prove the identity for one flavour N 2 = I first, by induction in N. For N = 1 
we have that 
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by adding v + 1 times the first column to the second column. Next we do the induction step, 

Lg(m) ... (l-r)- k Ll(m) ... (1 - t)-( n +V L" N+1 {m) 
Lg(xx) ... LUxi) ••• ^+i(xi) 
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(A.4) 
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L v N+l (x N ) 



Here we have expanded with respect to the last column and used the induction assumption for N, as 
well as the fact that the sub-determinant containing only x-variables of monic Laguerre polynomials 
can be replaced by the Vandermonde determinant. 

To get monic powers in the last column (except in the first element) we subsequently subtract 
multiples of columns from the left, using eq. (IX2|) . and we obtain 
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The first element in the last column now reads 
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As a last step we need to show that P(m) = t n+1 (1 — t)^^ n+1 ^ L^ +1 (™) . This relation holds due to 
the following identity [23], which can be easily proven by induction. It is expressed in terms of usual 
non-monic Laguerre polynomials 
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which finishes the first part of our proof. As a remark which is useful for the main text this identity 
is usually quoted as [23] 
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In the above it was not essential in the manipulation of columns that we had one mass flavour 
iV2 = 1 only. We can in fact do an inductive proof in the column number k for any N2 and N, 
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The start for k = 1 is trivially true in analogy to eq. (1A.3|) . The induction step from k to k + 1 easily 
follows by subtracting the left columns for I < k from column k + 1, and using again eq. (|A.7p for 
TV + 1 — ► fc + 1. Putting k = N + N 2 ends the proof. 

B An identity for Laguerre polynomials 

The relation we show here is given in eq. (|4.14p . 



where the right hand side is independent of k. It follows from a known identity eq. (4.4.1.14) in 

1 m 1 

Here (f3) m is the Pochhammer symbol. In choosing a = —j3 we obtain 

-L^x) = 2 _ if"*(x) . 



a! m v y £^ (m-*)!a! 

(_)m-j (a _ m + j)! Q _ m+ . 
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In the second step we have changed summation from i to j = m — i. Finally choosing a = m + k we 
obtain eq. (|B.1[) above. 
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